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Abstract 

This study delves into the analysis of the free vibration characteristics of a Timoshenko sandwich beam tapered 
along both width and thickness, subjected to an axially pulsing load on a Pasternak foundation. Integral to 
this investigation are the comprehensive considerations of energy expressions encompassing bending, axial 
deformation, kinetic energy, and strain energy induced by transverse shear stress. Subsequently, the 
differential equations governing the system, along with the associated boundary conditions, are deduced 
utilizing Hamilton's principle, followed by their non-dimensionalization. Leveraging series solutions from 
prior research that satisfy these equations and boundary conditions, each coordinate of the system is 
addressed. In the quest for a deeper understanding, Galerkin’s energy principle is employed to derive matrix 
expressions for key parameters such as mass and stiffness. Utilizing these matrices, Eigenvalues are computed 
to ascertain the natural frequencies of the system. The findings are presented through a series of graphical 
representations, elucidating the influence of various system parameters. 

Keywords: Timoshenko Sandwich Beam; Free Vibration; Natural Frequency; Taper Parameter; Temperature 
Gradient. 


1. Introduction 


In order to reduce noise and vibration, sandwich oncases where the core material exhibits viscoelastic 


beam are often used in the automobile, aerospace, and 
marine industries. The strength to weight ratio of 
Timoshenko sandwich beam is higher that of solid 
straight beams. In order to reduce noise and vibration, 
it is necessary to look at the needs in numerous 
industries and extensive study has been done on 
sandwich beams. Ray and Kar [1] investigated the 
instability of sandwich beams across various 
boundary conditions, scrutinizing parametric 
vibrations extensively. In a separate study, Kerwin et 
al. [2] delved into the damping effects on the free 
vibration of sandwich beams, particularly focusing 
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behavior. Meanwhile, Saito and Otomi [3] conducted 
an analysis on the stability of beams supported by 
viscoelastic supports, incorporating the influence of 
mass placement along the beam's length. 
Additionally, Kar and Sujata [4] observed a 
correlation between increasing taper parameters and 
decreased stability against periodic forces, as well as 
diminished static buckling loads. The core density 
parameter has no impact on the static but it does have 
a discernible effect on the beam’s sensitivity to 
pulsating force, either increasing or decreasing. 
Nayak and Dash [5] Observed that the static stability 
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of the beam increases as the pre-twist angle increases, 
but only for odd modes. Ansari et al. [6] studied a 
cantilever sandwich beam with tip mass and the effect 
of tip mass on the stability of the system. Lin and 
Chen [7] employed FEM to investigate the dynamic 
stability of a rotating sandwich beam. Their study 
showed that increasing the thickness of the 
viscoelastic layer and angular velocity led to 
improved dynamic stability. However, large setting 
angles were found to have a negative impact on 
dynamic stability. Nayak et al. [8] investigated a 
tapered beam of exponential variation, pre-twisted 
system on a Pasternak foundation supported 
viscoelastically and a temperature grade. Pradhan et 
al. [9] examined the stability of an uneven 
sandwiched beam on a foundation of Pasternak type 
with a temperature grade. The stability of a non- 
uniform uneven sandwich beam on a Pasternak base 
in a temperature environment was investigated by 
Pradhan et al. [10]. Malekzadeh et al [11] analysed 
on a functionally graded circular curved beam and 
found that the frequency parameter decreases as the 
opening angle increases. The dynamic Timoshenko 
beam subjected to axial periodic force was studied by 
Sabuncu and Evran [12]. The researchers concluded 
that as the distance between the centroid and the 
centre of flexure rises, the coupled bending dominant 
frequency grows. Ditaranto and Blasingame [13] 
found generalized result for composite loss factor and 
natural frequencies. The forced vibration was 
analysed by Mead and Marcus [14] for different 
boundary conditions by using the solutions provided 
by Ditaranto and observed different uncoupled 
complex modes. Banerjee [15] presented the 
problem-solving approach by dynamic stiffness 
method for sandwich beam. Nayak et al. [16] 
investigated a tapered beam of exponential variation, 
pre-twisted system on a Pasternak foundation 
supported viscoelastically and a temperature grade. A 
functionally graded, pre-twisted thick cantilever 
beam was the subject of an investigation by Mohanty 
et al. [17] into its vibration and dynamic stability. 
They discovered that a low index damages the 
system's stability when power levels rise. Liao & 
Huang [19] enriched the literature by exploring the 
parametric instability of spinning pre-twisted beams 
under periodic axial forces, contributing to a deeper 
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comprehension of the intricate phenomena dictating 
these structures. In a complementary study, Shiau & 
Tong [20] explored the intricate relationship between 
stability and response in rotating pre-twisted tapered 
blades, shedding light on the nuanced dynamics at 
play. Joubaneh et al. [21] conducted numerical and 
experimental analyses of free vibration in sandwich 
beams with tip masses using higher-order theory, 
highlighting the significant impact of tip masses on 
natural frequency. Freidani et al. [22] considered the 
consequence of a moving mass and curvature radius 
on dynamic deflection, observing — substantial 
increases with higher parameter values. As 
investigation of free vibration analysis of an 
asymmetric Timoshenko sandwich beam which has 
tapered along depth subjected to axial pulsating 
excitation is not available till date, the existing work 
shares out with the above said system configuration 
2. Problem Formulation 


Figure 1 shows a 3-layer Timoshenko tapered 
sandwiched beam with a periodic axial Joad. The 


periodic axial load is denoted by P(t)=R, +P cosat. 


: 
= 


\ 
4) 
Hii W | CROSS-SECTION A-A 
0.000.000 0.0. 0.0.0 0 Welrunsverse Deflection: Axial Deflection 


Figure 1 System Configuration 
2.1 Materials 


In this analysis, a three-layer sandwich beam is 
constructed using plywood as the viscoelastic middle 
layer, with structural steel and Cu-Al alloy being 
considered for the top and _ bottom sheets, 
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respectively. Using the right adhesive, the layers are 
joined. The materials used to build the beam are listed 
in Table 1 along with a variety of mechanical 
attributes. 


Table 1 Mechanical characteristics of the 
Considered Materials 
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beam. Its elastic modulus is far lower than those 
elastic layers. 


The total kinetic energy of the system can be written 
by 


DR ae . 42 
T = 5] cai’ + p,Au; + p,Aus+ p10) + pl,0;)dx (1) 
0 


The total potential energy of the system can be 
written by 


kapha + E10, + E,1,0,, 


yt} es 
29|(+K,AG, (Ww, -4) +K,A,G,y,' + K,A,G,(w, -6,) 


Materials | Density | highest Young’s 
(kg/m) | tensile modulus 
strength (MPa) 
(MPa) 
ASTM 7800 550 200*10° 
A36 Steel, 
plate 
Cu-Al 7850 470 170*10° 
Alloy 
CuAI8Fe3 
Materials | Density the Young’s 
(kg/m?) | highest | modulus (Pa) 
tensile 
strength 
(Pa) 
Plywood 650 13.8 620*10° 


The equations of motion for the beam are formulated 
under specific assumptions, which are as follows. 

1.The elastic layers of the skin adhere to Euler's 
theory, while the beam's core undergoes shear 
deformation theory 

2.The system's pre-twist angle is incredibly 
modest. 

3.It is presumed that there is no sliding at the 
intersection of the beam's several levels. 

4.It is possible to disregard the strain energy and 
extensional inertia of the core caused by in-plane 
stresses. 

5.The beam faces’ 
minimal. 

6.The deflection curves have short slopes. This is 
owing to the expectation that any deflection will 
be constrained by an increase in the stiffness of 
the system brought on by a high elastic modulus 
or a thicker viscoelastic layer. 

7.Jt is deemed possible to ignore the resulting force 
in the intermediate viscoelastic layer of the 


rotational inertia may be 
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1 to, Br , 
; aay a dx (2) 
The total work done of the system can be written by 


W = ; J P(t)w 2dx (3) 


Energy principle by Hamilton 6 J (T —V —W)dt =0, the 


following are the governing 


-Mw+ K,A,Gyw ,,. + K,A\G,4,,, + AGW an + K3A,G,03,, 
+G,bdw ,., = Bow POW yx 

(4) 
“(aA +a” p3Ay)ii, + (E,A + a7 Ey A3)Uy, ox 
K 
Ada uy - (1+ ahd, —(1+ a)hy0, |= 0 (3) 


2 


(-p,1,6, + E\LA, x. + Ki A\Gw , — Kj AGG, 


K,A,G, 2 = (6) 
era: A -(1+ a), +hyhy0; |=0 
CBee Bia AGN. — K;A,G3, ‘ 
— 52022 | 1° - (1+a)hyu, +hhyO, |= 0 7) 


and the boundary conditions for the system are 
-K, AG w, + K,A,GA, — K,A,G,w + 


8 
K,A,G,0, —G,bdw,, + P(t)w,, = 8) 
E,Au,, +a’ E,Au, , =0 (9) 
E,1,0,, =0 (10) 
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E,1,0,,, =0 


(11) 
The dimensionless form of total energy terms are as 
follows 


No BP 
M 


1 
[ AA, ae + 


10 10) 


lp —- 
P= 5 | ae 


a dx + (12) 


2 


hy BP ° (aA 
M 0 


10 


8) di 

é 0 
l+a?+ ~~ +h070, +h,0,0, + 2hh,0,0,0, 
—2Ni,9,9, 7 Nit, 2,9, ~ 2h,it,9,0, ~ 2ah,it,0,, 


~ 6, de+ 


K,A,G,(, 


1 
1 
ca 
3bdG L' 3 
SCE inti 
4Bh,, 
The total dimensionless work done becomes 


wh =| k@w'ae (13) 


10°10 0 


Te pes 
pee ue dx (14) 


The dimensionless governing equations for motion 
are 
—Mw + K,A\G\w + K,A\G,O, , + K3A,G3W ., + K3A,G304, 


+G, bdW — Bk(x)w- P(t Wx, =O 
(15) 
~(p\A, + + Ot” P3Ag id, + (BE, A; + a? By Ay), xy 
K,A ons —_ 
Aba +a’), —(+a)hO, - (1+ ands | =0 (16) 
Ah, 
ATO, + Eh 6, a t KjAGw, - KA, G 1, 
= RAO 7, [i 76, -(+ ayia, + ihs0 |= 0 (17) 
(-DyF,,, + ExT, + KyAsGyW , — KzA,G0, 
(18) 


me [i 0; (+ @)hyii, +h h,8, | = 


International Research Journal on Advanced Engineering Hub (IRJAEH) 


International Research Journal on Advanced Engineering Hub (IRJAEH) 


e ISSN: 2584-2137 
Vol. 02 Issue: 06 June 2024 
Page No: 1632- 1638 


https:/Arjaeh.com 
https://doi.org/10.47392/IRJAEH.2024.0224 


2.2 Solution by Variation Method 
The system's solution was estimated through the use 
of the following set of equations 


W(x, ft) {yi (19) 
i=l 


f 
uy (x,t) = >, GMix 


(20) 
k=e+1 
= & 
A, (x, t) a >, Imam (21) 
m=f +1 
_ h 
0; (x,t) = = Goo (22) 


o=gtl 


. The generalised coordinates are W, ,U,,, A», Aso 


The equations (19) - (22) provide the expressions for 
the shape functions that correspond to the pinned- 
pinned boundary condition as in [12]. 


w;(X) = sin(imx) (23) 
i, = cos(knx) (24) 
6, = sin(mrnx) (25) 
6; = sin(omx) (26) 


After non-dimensionalising, the extended Galerkin 
method is applied to the energy equations to obtain 
the following matrix equation 


[M , {O}+1[K]— FoLHI{O}— F, cos oF [H]{Q} ={9} 

(27) 
M, and K, are global mass matrix and global 
stiffness matrix correspondingly. 
RF, — fo 


28 
So at 


1635 


International Research Journal on Advanced Engineering Hub (IRJAEH) 
e ISSN: 2584-2137 

Vol. 02 Issue: 06 June 2024 

Page No: 1632- 1638 


https:/Arjaeh.com 
IRJAEH https://doi.org/10.47392/IRJAEH.2024.0224 


F, cos(@,t) (29) 2.3 Free Vibration Analysis 
For determination of the natural frequencies P(t) is 
occupied as zero (i.e., Fy =F, =0) in (27). The 


F(t) = 
10) 


The various sub matrices are : : ; 
natural frequencies can be obtained from the Eigen 


my = { Mw,w dx (30) Values of the mk matrix. Referring [13] from above 
0 the following equations can be derived 
3. Result and Discussion 


1 
Wasi = fA@Ai +00" p, Ag iy uy, dx (31) Numerical data are available for parameters such as 
: shear, thermal gradient, taper, spring characteristics, 
; and spring loss factors. Figures 2 to 8 depict the 
= APA, Oty, , 4% (32) influence of various system parameters on the 
0) 


33mn 
structure's free vibration. 
; 1.8 
Ma44op = J 273130505 )ae (3) 1. 
a‘ 1.2 
1 a ay ' ' a aS ' ' Bog 
[(A5K,A;,G,w; Wj ) + (A; K3A3,G3W, Wj )+ 
kui =| ae = (34) 0.6 
cal | 
0.3 
1 
K13mj = J As i GC, Fig ® (35) 1 02 03 04 05 068 O07 08 09 1 
d/l 
1 é = 
kisyy = [-AsKyAy,Gyw/ 0,4 (36) Figure 2 Response of @ for d/l 
10) 


Figure 3 depicts the impact of the g /j on the natural 


1 


k = [{A, (EA, + @°E,,A,, iy by }+ 
22k1 =| 


on (37) frequencies of the system. For larger vales of g /L > 
0{A,K,A,,G,0+ ay Uy, Uy, } \dx 


the natural frequencies of the system increase. 


1 
Keon, = [A+ @)A,K, A, Gh ptt, Oj, LE (38) 100 
0 
: SS 250) 
Kyat = [-G+ @)A,K,A,,Gh,,14,,0,,4% (39) 
2 200) 
1 rT ’ i 7, 7 
[(A,£,L9,,, 9, )+ AK, A,,G, A In) + 
Kam = f peer ; 5 1m 1 (40) 1B 150 
0 (AKA, Goh, AG An Ay lax 
1 100 
Ka10n = [46K A. GM, ,, O5,% (41) 50 
0 
1[(A, £41,0,, 0;,/) +(A,K;A;,,G,,,0. 
all ee = i Sere eee (a) 4 02 03 04 05 06 07 08 09 1 
0 (A,K,A,,.Gih, 0,5, )\dx B 
1 
Aug = J w/w dx (43) Figure 3 Response of @ for £ 
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Figure 3 depicts the impact of the width tapper pon 


the natural frequencies of the system. For larger vales 
of B° the natural frequencies of the system increase. 


83 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
G, /E, 


Figure 4 Response of @ for G;/E, 


Figure 4 The G,/E, outcome enhances the natural 


frequencies across all the three modes due to the 
improved rigidity of the system. 


4400 


120 


a 


$3 


Figure 6 Response of © for ¢3 


Figures 5 and 6 illustrate the change in natural 
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frequencies with increasing values of the thickness 

taper parameter. It's evident that as the taper 

parameter rises, the natural frequencies decrease. 

This decline is attributed to the reduction in cross- 

sectional area, which leads to a slenderizing effect on 

the beam and consequent deterioration in natural 
frequencies. 

Conclusion 

The findings derived from the comprehensive study 

shed light on the positive correlation between the 

width taper and the enhancement of the beam's 
natural frequency. Additionally, it was observed that 
the ratios involving the thickness of the shear layer to 
the length of the beam and the modulus of the shear 
layer to the modulus of elasticity of the elastic layer 
play pivotal roles in augmenting the natural 
frequencies. On the contrary, elevating the taper 

parameter along the thickness was found to have a 

detrimental effect on the natural frequency, while a 

corresponding increase along the width was 

associated with improvements in the system's natural 
frequencies. 
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